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A well known major difference between gravity and other field theories is that the former 
includes second order time derivatives in its usual Lagrangian formulation. In order to pass 
to a Hamiltonian formulation, Pirani et at circumvented this situation by subtracting 

a divergence term from the Lagrangian which includes second order time derivatives, leaving 
one with first order time derivatives of the fields in the Lagrangian and making it suitable 
for the Dirac treatment of constrained systems j^]. 1 More specifically, they split the d- 
dimensional Einstein-Hilbert (EH) action in the following way 2 

Seh = J d d x\] {-l) d - 1 g R = J d d xL rr - J d d xV a a , (1) 

where 



£rr = vH-irW^r^-r^) , 



V* = y/i-l^g^T^-g^j ■ 

and then neglect the contribution of the surface term V a a as its inclusion has no effect 
on equations of motion. As was emphasized in 1], a problem with this approach is that 
the part of the action remaining after elimination of the surface term, the rr-part, is not 
invariant under a general coordinate transformation. This can lead to inconsistencies in the 
Hamiltonian treatment of this system. 

The importance of a surface term for the Hamiltonian formulation of General Relativity 
based on the ADM slicing of spacetime was emphasized by Regge and Teitelboim 6] and in 
the path- integral approach by Hawking [3]. 

In two dimensions the gravitational action, given in (pQ), is a topological quantity. The 
equations of motion for g^ v in two dimensions are trivial identities, putting no restriction 
on the metric. No matter how the Cauchy problem is formulated, the gravitational fields 
are arbitrary functions of spacetime; as any possible configuration of metric extermizes the 
action, the latter is a constant. 

n 

However, in two dimensions not only does the Lrr part of ([1]) not vanish [8], but it also 
gives rise to a striking but consistent Hamiltonian treatment of two dimensional gravity. 



1 Another alternative way to circumvent this situation is by working with first order formalisms as briefly 
discussed later. 

2 is the Christoffel symbol, g = det Qu V , the signature is (+,—,—, ...), and the convention for defining 
the Riemann tensor is the one used in Q or Q]. 
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Applying the Dirac formalism to the Lrr leads to a gauge transformation [8] which is simply 
that one can add any arbitrary tensor to the metric tensor. This is consistent with the 
metric tensor being arbitrary. 

In contrast, the total two-dimensional Lagrangian of flTJ, when expressed in terms of 
tetrad variables, is a total derivative. This fact has eroneously led to the conclusion that a 
Hamiltonian formulation of two-dimensional gravity is impossible. 3 However, a Hamiltonian 
treatment of (PQ) is possible in the metric formulation on account of the inequivalence of the 
metric and tetrad formulations of gravity 9|. 4 

One approach to the Hamiltonian formulation of gravity which avoids the problem of 
second order time derivatives is to work with the first order formulations of gravity, such 
as the metric-affine connection formulation of Einstein or the tetrad-spin connection 
formulation ll| - 

For the two-dimensional case, however, the equivalence of the metric-affine connection or 
tetrad-spin coonection formulation with the original second order formulations is obscured, 
as the metric (tetrad) no longer uniquely determines the affine (spin) connection. This is 
explained in [l5( - 17|. Fortunately, there is still one more approach that is applicable 



to theories with higher-order time derivatives which involves reducing the order of time 

(see also (l9j). The 



derivatives by introducing extra fields, following Ostrogradsky 
application of such a formulation to two-dimensional gravity is the subject of this article. 

The Hamiltonian formulation of higher-order theories was considered more than one and a 
half centuries ago by Ostrogradsky [18J] for the case of non-singular systems, and was several 
times rediscovered by others. 5 Before generalization to singular systems, this ap proach was 
discussed for a few higher-order theories such as Podolsky electrodynamics 20|] - 22], and 



and was first applied to the EH action in 



scalar fields with higher-derivative cou plin gs 
four dimensions by Dutt and Dresden [24]. 

The first systematic generalization of the Ostrogradsky method to singular systems was 



3 This conclusion was supported by writing down the Lagrangian for metric gravity in special coordinate 
systems (such as the conformal frame), where the IT-part vanishes, and then wrongly generalizing this to 
any coordinate system. 

4 As the action is a total derivative in the tetrad formulation, any possible transformation of tetrads is a 
candidate for a gauge transformation. 

5 The criticism that can sometimes be found on the limitation of the original Ostrogradsky results to 
non-singular Lagrangians is not entirely correct because he clearly stated this restriction in his analysis. 
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considered by Gitman, Lyakhovich and Tyutin (GLT) |25| (see als o (26 
plied to the Hamiltonian formulation of generalized Yang-Mills theory 27 ] 
derivative gravity 29]. 



], and was ap- 



281 ] , and higher- 



By a suitable introduction of extra fields, this generalization of the Ostrogradsky method 
allows one to reformulate a problem with higher derivatives in such a way that the Dirac pro- 
cedure which was originally capable of handling only theories with first-order time derivatives 
[2I can be used for singular, in particular, gauge theories with higher order time derivatives. 
Details of this generalization and its variations appearing in the literature 30(], 3jj, and 32] 
are beyond the scope of this article. We instead provide an illustration of how it can be used 
by considering the EH action in d dimensions and then specializing to d = 2 dimensions. 

The Lagrangian of the (metric) EH action depends on second derivatives of the fields 



Lg — L g {g a /3] 9a/3,0] 9a/3,k] 9a/3,00] 9a/3,0k] 9a/3,km) 



(2) 



where we have separated time derivatives of the fields to make subsequent discussion more 
transparent. 

If, following GLT, the additional variables 



Gaf3 — g a /3 ; V a/ 3 — G 



ai3 



(3) 



are introduced, the Lagrangian ([2]) can be represented in the following way 



L V q — L V „ {g a (3\ G a /3] g a f3,k'i V a /3] Ga(3,k\ ga/3,km) 



(4) 



where only first-order time derivative of fields appears. As it will be shown below, this 
Lagrangian describes the same dynamical system as that of ([2]) if it is supplemented by the 
conditions of (ED 



S 



L v g d d x 



L v g + tt^ (g aP - G a/3 ) + IT* (G. 



a/3 — V a/ 3 



d d x, 



(5) 



The Lagrange multipliers ir a ^ and U a/3 act as momenta conjugate to g a p and G a p respec- 
tively. Variation of this action results in the following equations of motion 



5L v a dVL dL v a dL v a , 

7 = o 9 k - — — + d k d m — S- 7r = , (8) 

v9ap O9o.p 09af3,k C>9oif3,km 

$K dL v dU a R . R 

9 — = i o k - 7c al3 - n Q/3 = o , (9) 

5G a/ 3 dG a/ 3 dG a/ 3 t k 

$L v n dL v n R 

-^ = 7-^-n Q/3 = 0, (10) 

bv a/ 3 dv a/ 3 

provided all variations vanish on the boundary. These are equivalent to the Lagrange equa- 
tions following from variation of the original EH action (T5]) 6 

5L„ dL„ „ dL„ „ dL G 



9 - 9 d n ^- - a ^ 9 



^9a/3 dg a p dg af 3fi dg a p^ 
dod OJ ^- + d d kJ ^- + d k d m -^- = . (11) 

CflW?,00 OgaPfik Ogap^km 

This is easy to verify by differentiating ([TO]) with respect to time and substituting I1 Q/3 into 
(JSJ). We differentiate again with respect to time to get 

dL v n dVi dL v n 

* = ^—^ ~ d Q d k —^- - d d -^ . (12) 

OLr a p OLr a p^ k OVa/3 

Substituting of n alS into (jHJ) and supplementing the result with solutions of (jEJ) and (J7j) 
finally establishes equivalence. The advantage of the new Lagrangian is that it contains 
only first-order time derivative of fields, permitting us to write the Hamiltonian in the usual 

way 

H = ir aP g aP + Ii aP G a(3 - L v g , (13) 

leading to 

H = T^G aP + U^v a0 -L v g . (14) 

This can be verified by a simple rearrangement of the terms in V. The fundamental Poisson 
brackets (PB) associated with this Hamiltonian are 

{g^ iT) = {G a ,, Tin = i (6*% + 6$%) , (15) 

and the PB for any two functionals A and B of the canonical variables is 

r , 8 A 5B 8 A 5B 

{A, B} = - — — + ^ — — {A <— » B) . 16 



6 Note that there is no coefficient "2" in the fifth term of (fTTj) since the symmetry of double derivatives 
dod k has been implicitly taken into account in @ (for an explicit form of see (l2l)|) below). 
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As an indication of the singular nature of this Lagrangian, we observe that the equations of 
motion following from the Hamiltonian formulation are 

5H 

9af3 = {gap, H} = ^ = G a p , (17) 

5H 

G a( s = {G a/ 3, H} = -^-^ = v a/ 3 , (18) 
^ = {^ iH } = ~ = p-, (19) 

il^ = {n^,H} = -£L = ^L, (20) 

which are not equivalent to the equations of motion following from the Lagrangian formu- 
lation. This is because equation ffTUl) 

dL v 

- U a/3 = (21) 

ov a p 

is missing. If the Lagrangian were non-singular (Hessian were not zero), one would be 
able to solve (I2T!) for v a p and substitute the solution back into ({TBI to get a consistent 
Hamiltonian formulation for the dynamics of the two pairs of conjugate variables (g Q/ 3, 7r a/3 ^ 
and (G a p,Tl a ^ , without any of constraints. However, for the EH Lagrangian, because it is 
linear in the second-order derivatives, Eq. (I2TI) cannot be solved. Therefore Eq. ([21!) has to 
be supplemented to the Hamiltonian formulation as a set of primary constraints. 

We now demonstrate this explicitly, starting in d dimensions, and then switching to two 
dimensions, which is the main concern of this article. If the Riemann tensor is written 



explicitly in terms of the metric tensor and its derivatives 24fl . the EH Lagrangian splits 
into two parts in the following way 



where 



and 



-iy-ig R = A^g + 5*^, 7 v,p , (22) 



= sj{-iy-ig (g a ^ - g a(i g^) (23) 



+2g ap g f3u g^ + Ag^g^g^ - 4# a V V*J • (24) 
After introducing the additional variables of ([3]), the Lagrangian becomes 

L = T^g aP + Ji^G a p-H, (25) 

6 



with 

H = n at3 G a(3 + v af3 (iV* - A at300 ) - (A a ? ok + A at3k0 ) G af3 , k - A^ km g af3 , km 
- B a ^ 00 G a pG^ - (B a ^ ok + B ^ k0 ) G aP g^ k - B a ^ km g a ^ k g^ m . (26) 



Here it can be explicitly seen how the second term in this expression leads to the aforemen- 
tioned primary constraints of (TTD1) 



where A al3 ^ u is given by fl23|) . At this point we set d = 2. By making use of the fact that in 
two dimensions 

one obtains the following constraints from (J2j 



p 11 = n 11 - -L= « o , p 00 = n 00 « o , p 01 = n 01 « o . (28) 



-9 

These primary constraints P a @ obviously constitute a first class system with a simple algebra 
of PBs 

{P Q/3 ,P^} = 0. (29) 

The secondary constraints are determined by requiring that the primary constraints be 
preserved in time 

pa/3 = j po H\=S«P = Q) (30) 

leading to the following secondary constraints 

1 a 00 1 a 01 a 00 

S 00 = -vr 00 - i-?=G u + £-7=011,1 + -7=0)1,1 , (31) 



2V=9 2^ ' V=0 

1 O 01 1 fl 00 



2^=9 2^=0' 

1 O 11 1 fl 01 

5 11 = -vr 11 - ^-7=^11 - ^-7=^1 • (33) 



2V=9 2^' 

It is not difficult to show that these new constraints ( 13TH331) commute with the set of primary 
constraints 

{p af3 ,S^} = 0, (34) 
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and the PBs among secondary constraints S a/3 also vanish 

{S a(3 ,S^} = 0. 



(35) 



Thus all primary and secondary constraints are first class. Using (13111331) one can write the 
Hamiltonian (1261) in the following way 



with H including only spatial derivatives 



H = v a(3 P a P - G aB S al3 + H 



a(3>- 



(36) 



H 



A a/311 r>af3(ivll ri 
~ A 9a(3,U D g a /3,l9iiv,l 

( fl'oo.i — 



(A * 01 + A * 10 ) G 



a/3 



01 



(37) 



V V-9 /,i 

The next step is to consider the time derivative of the secondary constraints to see if any 
new constraints arise 

S a ? = {S a ?, H) . (38) 

Using the form of the Hamiltonian in fl36l) and the fact that primary and secondary con- 
straints are all first class, these PBs are equal to jvr Q/3 ,if j. These latter brackets can be 
obtained using f|T3|) and assuming that the fields and their spatial derivatives vanish rapidly 
at infinity. An alternative is to treat if as a spatial surface term and totally neglect it. Either 
method shows that these brackets vanish, leading to the closure of the Dirac procedure with 
six first class constraints with zero PBs among them. As there are six independent fields g a p 
and G a p, and there are six first class constraints, there is zero degrees of freedom. Using the 



Castellani procedure 



331 ]. one can restore the gauge transformation of fields by building the 



gauge generator G (e) corresponding to the complete set of first class constraints. In much 
the same way as in 



341 ]. |35|, and [17]], we find that the gauge generator G (e) is 
G(e) = J dx (- £aP [P aP , H) + e aP P a ^ 



(39) 



This gives the following gauge transformations with arbitrary gauge parameters e a f : 



^9 a p = {g a p,G(e)} = e 



a/3 



and 



8G a/ 3 — {G a/ 3, G(e)} — Saf, 



(40) 



(41) 
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This is consistent with (|3]). The transformations of (|40j) are exactly the same as the trans- 
formations obtained in 8], where only the Lpr part of the EH Lagrangian was considered. 
However, there is a difference in the number of fields and the structure of the constraints in 
these two methods. 7 

We now investigate the application of the Ostrogradsky method to tetrad gravity. This 
merits interest since, as has already been discussed, the Lagrangian for tetrad gravity is a 
pure surface term in two dimensions, making its treatment using the method employed by 
Pirani et al impossible. 

If one substitutes into the EH Lagrangian ( 1221) the expression for in terms of tetrads 

= effiriab , (42) 

one obtains the second order Lagrangian for gravity in terms of tetrads, which can be 
analyzed using the Ostrogradsky method. 8 The result of this substitution can be written in 



the compact form 36] 



L e = (2e ab e^ a e bp>v ) ^ = 2e ah e^ (e£e ta - eje^e^) , (43) 

where e is the totally antisymmetric tensor e m = = l. 9 Introducing, as in the metric 

case, additional variables 

e; = v ; = e; (44) 

and performing manipulations similar to those used in obtaining ffl3|) . we get 

L: = ^ + W a E^-H } (45) 

with 

H = tfE* + ~ K- (46) 
We define the following fundamental PBs 

{e a ,,<} = {E;,Tl»} = S a b 6; (47) 



7 In passing, we note that there are a few possible variants of the generalized Ostrogradsky method that 
lead to the same transformations. For example, it is enough to introduce extra fields only for those 
g a f3 components that have second-order time derivatives (only g\\ in two dimensions). In this case, the 
analysis leads to four first class constraints (three primary and one secondary) for the four fields which 
again leads to zero degrees of freedom and the same gauge tramsformation as fTO)) . 

8 rjab - Minkowsky metric, [i = 0, d — 1 are world indices and a = 0, ...d — 1 are tetrad indices. 

9 () brackets distinguish explicit values of tetrad indices. 
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and have 



LI = 2e ab e° a v bl + 2e ab e 1 a E blA - 2e ab e° a E b0A - 2e ab e 1 a e b0 , n 



- 2e ab ele° c E c a E bl - 2e al 'e^,^ + 2e ab e:e° c E c a e m + 2e ab ele x ~ c 



(48) 



(49) 



After rearrangement of the terms in (1481) . one obtains 

H = W a v; - 2e ab e° a v bl + + 2e ab ey c E c a E bl - 2e ab eZe°e m>1 

+ 2e ab ey c e^E m - [2e ab (e l a e b ^ + e\E bl - e^o)^ • 
The first two terms in this expression give rise to a set of primary constraints P% 

P(o) = n° ) ) — n°!) , P( ) = II( ) + 2e°!) , P ( \) = n^) + 2e° 0) . (50) 

The PBs among these four primary constraints P£ are obviously zero. Conservation of 
primary constraints in time leads to the secondary constraints S ( 



a 



P^ = {P^,H} = S^ = 0. 



(51) 



They are given by 



^(0) - _7r (0) + 2e (0) e (0)^l + ^(l^O)"^ 

,9 1 (0) , p i (1) , o / _ N (1) 
+ ze (i) e (o) e i,i + Ze (l) e (l) e l,l + Z \^ e (l) e (l) e (o) e (o) J e o,i 



Q° — -rr _J_ 9«° «° /T^ 1 ) _L 9«° «° f?(°) 



,9 1 (0) , 9 1 (1) , 9 / _ \ (0) 
+ Ze (0) e (0) e l,l + Ze (0) e (l) e l,l + Z V (°) (°) (!) (!) / O- 1 ' 

Cl _ _ 7r l , O p l „0 R (1) , 9 1 „0 R (0)_ 9 1 „0 (1)_ 1 (0) 
D a — + Ze (0) e a- C/ 1 + Ze (l) e a- C/ 1 /e (0) e a e 0,l /e (l) e a e 0,l 



It is fairly easy to demonstrate that 



(52) 



(53) 



for any pair of primary and secondary constraints, but calculation of PBs among secondary 
constraints is slightly more involved, and requires the use of test functions, leading to 



{S£, Sfi = . 



(54) 



The Hamiltonian ( 1491) can be expressed in terms of a linear combination of constraints plus 
a spatial derivative term 



H = v«P» - E*S£ - 



2e ab (ele m + e\E bl - e° a E, 
10 



60 



(55) 



Again, there are no tertiary constraints because S% = {S£, H} = due to eqs. (155]) and 
(l54j) . The Dirac procedure is closed, with ei ght first class constraints leading to zero degrees 



of freedom. Using the Castellani procedure [33], the generator of gauge transformation is 



Ge (e) = Jdx {-el {P^ H} + e^) , (56) 
where e a ^ are gauge parameters. This generator leads to the gauge transformation 

6 e % = {e»,G e (e)} = e* i (57) 

for the tetrad fields. 

As discussed in the introduction, the transformation (|571) . as well as the transformation 
found for the metric formulation of the EH action ( HUl) . is consistent with triviality of the 
equations of motion for the EH action in two dimensions. The number of gauge parameters 
for the metric case is three, and for the tetrad case four, the same as the number of inde- 
pendent fields in each case, leading to zero number of degrees of freedom according to the 
standard way of counting the degrees of freedom. 

As has been discussed elsewhere, the tetrad representation of a spacetime leads uniquely 
to the metric representation [9]. However, the converse is not always true as it is not possible 
to uniquely determine a tetrad representation from a metric representation. For the same 
reason, a gauge transformation in the tetrad formulation corresponds to a unique gauge 
transformation in the metric formulation 

e^v = S e g^ = 5 e (riabe^ety = ^e^e^ + r] ab e h v 8 e e a lx = e bfl e b u + e hv e\ , (58) 

while a gauge transformation in the metric formulation does not lead to a unique gauge 
transformation in the tetrad formalism even if a unique tetrad system is specified. This is 
because the three equations of fl58|) can not be solved for the four unknowns e a v in terms 
of e M „. In 2D the Lagrangian for tetrad gravity is a pure surface term but its Hamiltonian 
treatment using GLT generalization is possible. This can provide an alternative approach to 
demonstrate unequivalence of 3D tetrad gravity and the Chern-Simons theory [37| (which 
differ by a total derivative) based on their Hamiltonian formulations. 

To conclude, the Hamiltonian formulation of the two-dimensional EH action as a higher- 
derivative theory leads to a consistent structure of constraints, and a vanishing number of 
degrees of freedom. The gauge transformations of fHUj) and floTj) are different from linearized 
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coordinate transformations. 10 The number of constraints required to reproduce linearized co- 
ordinate transformations as gauge transformation of the two-dimensional EH action, results 



in there being a negative number of degrees of freedom 38] which is clearly unacceptable. 
Obtaining diffeomorphism invariance as a gauge transformation leads to discrepancies in 
the number of degrees of freedom and the number of first class constraints appearing in 
the Hamiltonian analysis of other two-dimensional models. As an example, a Hamiltonian 
analysis of a scalar field in curved spacetime —gg^d^d^ gives five first class constraints 



for the four fields ((^,0) as reported in [39|, leading to minus one number of degrees of 



freedom if the diffeomorphism and Weyl invariances are to be gauge symmetries. However, 
as will be reported elsewhere, the treatment of this model using the Dirac procedure removes 
the contradiction arising from having a negative number of degrees of freedom and leads to 
a gauge transformation distinct from the linearized coordinate transformation. 

In this paper the method of GLT was emlpoyed in which only the order of temporal 
derivarives was decreased (the same as was done by Dutt and Dresden). An alternative 
way is to introduce extra fields also for spatial derivatives of fields, an approach that may 
be called a covariant Ostrogradsky method. Recently, the lowering of order of derivatives 
in covariant form was used for three dimensional topologically massive electrodynamics 
(TME) [40] to construct its first order formulation. This covariant Ostrogradsky method 



leads to consistent results, despite of the conclusion of 



and, in particular, preserves the 



gauge invariance which was explicitly demonstrated using the Dirac formalism in 42]. 
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